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ABSTRACT 

The similarity of quasar line spectra has been taken as an indication that the 
emission line clouds have preferred parameters, suggesting that the environment 
is subject to a fine tuning process. We show here that the observed spectrum is 
a natural consequence of powerful selection effects. We computed a large grid of 
photoionization models covering the widest possible range of cloud gas density 
and distance from the central continuum source. For each line only a narrow 
range of density and distance from the continuum source results in maximum 
reprocessing efficiency, corresponding to "locally optimally emitting clouds" 
(LOG). These parameters depend on the ionization and excitation potentials 
of the line, and its thermalization density. The mean QSO line spectrum can 
be reproduced by simply adding together the full family of clouds, with an 
appropriate covering fraction distribution. The observed quasar spectrum is a 
natural consequence of the ability of various clouds to reprocess the underlying 
continuum, and can arise in a chaotic environment with no preferred pressure, 
gas density, or ionization parameter. 



Subject headings: quasars: emission lines 
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1. Introduction 

The question of the origin of quasar emission hne regions is central to cahbrating these 
objects as probes of the Universe at the epoch of galaxy formation. Once we understand 
the spectrum in detail there is the potential of measuring their distances (and hence Qo 
and Ho) and of determining the composition of the gas, thus constraining nucleosynthesis 
during galaxy formation. 

Although the promise is great, progress has been slow. The lines are produced by a 
dilute gas far from LTE. Its spectrum is sensitive to line formation details, particularly the 
cloud gas density and flux of radiation (Netzer 1992). These are often combined into a 
single "ionization parameter" , the ratio of photon to gas densities. The earhest observations 
showed that the family of quasars have surprisingly similar emission line spectra (Davidson 
&; Netzer 1979), and this is often taken as evidence that a process adjusts clouds to a 
certain ionization parameter (Davidson 1977, Krolik, McKee, & Tarter 1981). The possible 
presence of an unspecified tuning mechanism made it impossible to know how to interpret 
correlations between emission line parameters and other properties {e.g., Baldwin, Wampler, 
& Gaskell 1989). Finely tuned cloud parameters suggest a highly ordered environment, not 
the chaotic one that might naturally result during the initial stages of galaxy formation. 

More recent work by Rees, Netzer, & Ferland (1989), Krolik et al. (1991), Goad, 
O'Brien, & Gondhalekar (1993), and O'Brien, Goad, & Gondhalekar (1994, 1995) have 
investigated distributed emission line cloud properties, but still assumed the existence of a 
single cloud pressure law with radius. Here we use a much more general approach to show 
that a powerful selection effect operates, so that a random mix of cloud properties will 
produce the spectrum of a typical QSO. There is no spectroscopic evidence for finely tuned 
cloud parameters. The spectrum of a typical quasar is consistent with spectral formation in 
an environment with no preferred pressure, gas density, or ionization parameter. 



2. A Single Cloud 

We model individual clouds using the following assumptions: a single column density of 
N(H) = 10^^ cm~^, constant hydrogen density throughout each cloud, and solar metallicity. 
Although below we show results for a single cloud column density, we expect there to be 
a broad range of column densities present in the BLR. Preliminary investigations into the 
sensitivity of the results to the presence of a range of column densities indicate that they 
are not strongly affected as long as optically thin clouds do not dominate the total line 
emission. With the exception of N V A1240, most of the prominent emission lines are not 
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strongly affected by our choice of metallicity (see Hamann & Ferland 1993). 

The ionizing continuum shape chosen was a combination of a oc v~^-^exp{—hi'/kTcut) 
UV-bump with an X-ray power law of the form f,^ oc v"^ spanning 13.6 eV to 100 keV. 
The UV-bump cutoff temperature, Tcuti was chosen such that the energy in the UV-bump 
peaked at 48 eV. The two continuum components were combined using a typical QSO UV 
to X-ray logrithmic spectral slope of a^x = —1.4 (note the explicit minus sign). Korista et 
al. (1996) will explore in detail the dependencies of the emitted line spectrum on continuum 
shape, chemical abundances, and cloud column density. 

The remaining parameters are the cloud particle density n{H) (cm~^) and the photon 
flux of H-ionizing radiation $(i?) (cm~^ s~^). For an isotropic continuum source ^{H) 
is a direct measure of the distance between the clouds and the continuum source r. The 
variables r and n{H) directly describe the overall structure of the BLR, and we do not 
combine them into an ionization parameter U{H) = ^{H)/n{H)c. 

3. Locally Optimally-emitting Clouds (LOG) 

Figure 1 shows the results of a large series of photoionization calculations in which 
$(i/) and n{H) were varied. Lines are presented as equivalent widths for full geometrical 
coverage, referred to the incident continuum at 1215 A. This equivalent width is a direct 
measure of the cloud's reprocessing efficiency. 

CoUisionally excited lines such as C IV A1549 show a band of efficient reprocessing 
running at constant U{H) along a diagonal ridge from high ^{H) and n{H), to low ^{H) 
and n{H). Along this ridge, corresponding to log U{H) ^ —1.5 for A1549, the gentle 
decrease in the Wx at the high densities is caused by line thermalization — high density 
clouds are continuum sources. The line equivalent width decreases sharply when moving 
orthogonal to the ridge because U (H) is either too low (lower right) or too high (upper left) 
to produce the line efficiently. Clouds with log U{H) ^ 0.2 are optically thin so reprocess 
little of the incident continuum. C IV can be contrasted with a low ionization line such 
as Mg II A2798, whose peak Wx is shifted to lower U{H), or to the high ionization O VI 
A1034, shifted to higher U{H) 1). Note that the contours in Figure 1 are logarithmic; in 
linear space they are sharply peaked. We then see that the line formation radius naturally 
depends on the ionization potential (IP), modulated by thermalization density; low IP 
lines form at large radii and high IP at small radii, a result consistent with line-continuum 
reverberation measurements (Peterson 1993). 

Recombination lines of H° and He"*" emit over a wider area on the density-flux 
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plane, including the low ^{H) - high n{H) regions since both ions still exist under these 
conditions. On the whole, Lycu forms closer in than B.f3. This is because Lya is a resonance 
line and is collisionally suppressed at high densities, while B.(3 is suppressed at high flux due 
to larger excited state optical depths (Ferland, Netzer, & Shields 1979). This will make Lya 
respond first to continuum variations and make its line profile broader in the wings if cloud 
motions are purely virial. The first effect is observed (Clavel et al. 1991; Peterson et al. 
1991); the second effect is often but not always observed (Zheng 1992; Netzer et al. 1995). 
For He II A1640, the region of most efficient reprocessing is a diagonal plateau (rather than 
a ridge), peaking and saturating at high density and flux, but also diminishing at very low 
ionization parameter. This peak and much of the plateau are at signiflcantly larger flux 
than for H/?, thus A1640 is formed in gas lying closer to the continuum source. This is also 
consistent with the reverberation measurements (Korista et al. 1995). 

4. The Integrated LOC Spectrum 

The equivalent widths for full coverage, as plotted in Figure 1, represent the efficiency 
with which the individual LOC clouds can reprocess ionizing continuum radiation into line 
radiation. These reprocessing efficiencies act as a filter whose output is the QSO's emission 
line spectrum. Each line can be formed efficiently only at a particular location in the 
density-fiux plane. This location is different for different lines. If the individual clouds are 
distributed over a wide range on the density-ffux plane, then most emission lines will be 
produced with high efficiency, so long as that component has a significant covering fraction. 

The simplest description of this type of system is to take, for each separate line, the 
extreme case of considering only emission from a cloud with the optimal density and fiux 
parameters for the line in question. Under these conditions the relative intensity of each 
emission line would be directly proportional to the maximum reprocessing efficiency, for 
uniform covering fraction. Column 2 in Table 1 presents a mean observed QSO spectrum, 
and in column 3 are the results of this simple model. The two are in good agreement. 

A more realistic refinement is to assume some distribution of clouds on the density-fiux 
plane and use that as a weighting function to integrate over the distribution of reprocessing 
efficiencies. The limits of this integration are set by physical or observational considerations: 
clouds at large distances from the continuum source (log$(i?) <18) will form graphite 
grains (Sanders et al. 1989) and therefore have very low emissivity (Netzer & Laor 1993); 
very low density clouds {n{H) < 10^ cm~^) are ruled out by the absence of broad forbidden 
lines, such as [O III] A4363. Clouds at very small radii or with very high densities are either 
at the Compton temperature or are so dense that they are continuum rather than emission 
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line sources. Thus to produce the emission hne spectrum, one simply integrates the line 
emission from clouds along each axis in Figure 1 {logn{H) > 8, log$(i?) > 18). 

Remaining to be specified are the distribution of clouds as functions of distance from 
the continuum source and gas density, appearing in the definition of the emission line 
luminosity. 



where F(r) is the emission line fiux of a single cloud at radius r, and /(r) and g{n) are 
the cloud covering fractions with radius and gas density, respectively. Column 4 of Table 1 
shows the results for the simple case of /(r) oc constant and g{n) oc n~^; the results are 
sensitive to these distributions and will be discussed in Korista et al. (1996). The agreement 
with the observations again is quite good, with the exception of N V A1240 (underpredicted, 
as usual for solar abundances, see Hamann and Ferland 1993). Some of the other weaker 
lines in Table 1 whose predicted relative intensities lie near the bottom of the observed 
distribution should also be larger in gas of higher metallicity expected in QSOs {e.g., Si III], 
Si IV, Al III). The Lja/B.(3 ratio is still problematical, with large uncertainties in intrinsic 
reddening and radiative transfer in the Balmer lines (Netzer et al. 1995). In this example, a 
total geometrical covering fraction of 18% will produce the typical l^A(LyQ;) = 100 A. 



These results show that line emission from the BLR is dominated by powerful selection 
effects. Individual BLR clouds can be thought of as just machines for reprocessing radiation. 
As long as there are enough clouds at the correct radius and with the correct gas density 
to efficiently form a given line, the line will be formed with a relative strength which turns 
out to be very similar to the one actually observed. Thus there is no spectroscopic evidence 
for preferred cloud parameters. Observed QSO spectra carry information mostly about 
global properties such as the continuum shape, the clouds' overall metallicity, and their 
distribution, rather than about the details of individual clouds. 

This in turn shows that a jumbled or chaotic cloud environment can be the source of 
the lines. The spectrum of a typical quasar may be produced in a scrambled environment 
similar to that which might occur in the center of a proto-galaxy, with rapid star formation, 
evolution, and supernova activity. 

A ramification of the fact that lines form at different places is that we can qualitatively 
reproduce many of the observed variability and line width differences. This picture also 
provides a natural explanation for why the sizes of the emission line regions scale with 
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luminosity (Peterson 1993); lines naturally form at the appropriate ionizing fluxes, and the 
corresponding radii scale with luminosity. 

Many aspects of the LOG picture are under investigation: sensitivity of results to 
ionizing continuum shape, metallicity, cloud covering factor distributions with radius and 
gas density, resulting line profiles and reverberation. The likely presence of a broad range 
of cloud column densities converts the LOG plane to a cube {^{H),n{H), N{H)). 

If the LOG picture proves to be correct then our understanding of the BLR will have 
undergone the same change as did our view of the Narrow Line Region when it was realized 
that all of the forbidden lines were radiating near their critical densities (Filippenko & 
Halpern 1984). Simple averaging, not a hidden hand or unknown physics, is at work. 

We are grateful to Mark Bottorff, Moshe Elitzur, and the referee Simon Morris 
for helpful comments. This work was supported by the NSF (AST 93-19034), NASA 
(NAGW-3315), and STScI grant GO-2306. 



TABLE 1 

Observed and Predicted Line Intensities 



Emission Line 
(1) 


Observed Intensity" 
(2) 


Maximum Reprocessing 

(3) 


LOG Integration'' 
(4) 


VI A1034+Ly/?A1026 


0.1-0.3 


0.28 


0.16 


Lya A1216 


1.00 


1.00 


1.00 


N V A1240 


0.1-0.3 


0.06 


0.04 


Si IV A1397+0 IV] A1402 


0.08-0.24 


0.08 


0.06 


G IV A1549 


0.4-0.6 


0.54 


0.57 


He II A1640 + III] A1666 


0.09-0.2 


0.11 


0.14 


G IIl]+Si III]+A1 III A1900 


0.15-0.3 


0.28 


0.12 


Mg II A2798 


0.15-0.3 


0.38 


0.34 


H/3 A4861 


0.07-0.2 


0.08 


0.09 



"Intensity relative to Lya A1216, combining data from Baldwin, Wampler, & Gaskell (1989), Boyle 
(1990), Cristiani & Vio (1990), Francis et al. (1991), Laor et al. 1995, Netzer et al. (1995), and Weymann 
et al. (1991). 

^Co-addition of emission from clouds as described in the text. 
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Fig. 1. — Cloud reprocessing efficiency is shown as a function of the hydrogen density and 
flux of ionizing photons. Line strengths are expressed as equivalent widths referenced to the 
incident continuum at 1215 A. Each bold line is 1 dex, and dotted lines represent 0.2 dex 
steps. The smallest (1 A) and largest decade contours are labeled for each hne; note that 
the contour values continue to increase beyond the largest labeled decade in all cases but H/3 
and He ll. The star is a reference point marking the "standard BLR" parameters discussed 
by Davidson and Netzer (1979). 



